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Abstract 



We have calculated Heisenberg exchange parameters for bcc-Fe, fcc-Co, 
and fcc-Ni using the non-relativistic spin-polarized Green function technique 
within the tight-binding linear muffin-tin orbital method and by employing 
the magnetic force theorem to calculate total energy changes associated with 
a local rotation of magnetization directions. We have also determined spin- 
wave stiffness constants and found the dispersion curves for metals in question 
employing the Fourier transform of calculated Heisenberg exchange param- 
eters. Detailed analysis of convergence properties of the underlying lattice 
sums was carried out and a regularization procedure for calculation of the 
spin-wave stiffness constant was suggested. Curie temperatures were calcu- 
lated both in the mean-field approximation and within the Green function 
random phase approximation. The latter results were found to be in a better 
agreement with available experimental data. 
PACS numbers: 71.15.-m, 75.10.-b, 75.30.Ds 

I. INTRODUCTION 

The quantitative description of thermodynamic properties of magnetic metals is chal- 
lenging solid state theorists since decades. Thanks to the development of density functional 
theory and its implementation into ab initio computational schemes, an excellent under- 
standing of their ground state (i.e., at T = K) has been achieved. On the other hand, 
most of the progress towards a description of magnetic metals at non-zero temperature has 
been based upon models in which the electronic structure is oversimplified and described in 
terms of empirical parameters. Although this approach has the great merit of emphasizing 
the relevant mechanisms and concepts, it cannot properly take into account the complex 
details of the electronic structure and is therefore unable to yield quantitative predictions 
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of the relevant physical quantities such as spin-wave stiffness, Curie temperature Tc, etc., 
for comparison with experimental data. 

It is therefore of a great importance to develop an ab initio, parameter-free, scheme for 
the description of ferromagnetic metals at T > K. Such an approach must be able to go 
beyond the ground state and to take into account excited states, in particular the magnetic 
excitations responsible for the decrease of the magnetization with temperature and for the 
phase transition at T = T^. Although density functional theory can be formally extended to 
non-zero temperature, there exists at present no practical scheme allowing to implement it. 
One therefore has to rely on approximate approaches. The approximations to be performed 
must be chosen on the basis of physical arguments. 

In itinerant ferromagnets, it is well known that magnetic excitations are basically of 
two different types: (i) Stoner excitations, in which an electron is excited from an occupied 
state of the majority-spin band to an empty state of the minority-spin band and creates 
an electron-hole pair of triplet spin. They are associated with longitudinal fluctuations of 
the magnetization; (ii) the spin-waves or magnons, which correspond to collective trans- 
verse fluctuations of the direction of the magnetization. Near the bottom of the excitation 
spectrum, the density of states of magnons is considerably larger than that of corresponding 
Stoner excitations, so that the thermodynamics in the low-temperature regime is completely 
dominated by magnons and Stoner excitations can be neglected. Therefore it seems reason- 
able to extend this approximation up to the Curie temperature, and to estimate the latter 
by neglecting Stoner excitations. This is a good approximation for ferromagnets with a 
large exchange splitting such as Fe and Co, but it is less justified for Ni which has a small 
exchange splitting. 

The purpose of the present paper is to describe the spin-wave properties of transition 
metal itinerant ferromagnets at ab initio level. With thermodynamic properties in mind, 
we are primarily interested in the long-wavelength magnons with the lowest energy. We 
shall adopt the adiabatic approximation in which the precession of the magnetization due 
to a spin-wave is neglected when calculating the associated change of electronic energy. 
Clearly, the condition of validity of this approximation is that the precession time of the 
magnetization should be large as compared to characteristic times of electronic motion, 
namely, the hopping time of an electron from a given site to a neighboring one, and the 
precession time of the spin of an electron subject to the exchange field. In other words, 
the spin-wave energies should be small as compared to the band width and to the exchange 
splitting. This approximation becomes exact in the limit of long wavelength magnons, so 
that the spin- wave stiffnesses constant calculated in this way are in principle exact. 

This procedure corresponds to a mapping of the itinerant electron system onto an effec- 
tive Heisenberg Hamiltonian with classical spins 

HeS = — X! "^ij^i ' ' (1) 

where Jij is the exchange interaction energy between two particular sites (i, j), and ej, e^- are 
unit vectors pointing in the direction of local magnetic moments at sites (^, j), respectively. 
The same point of view has been adopted in various papers recently published on the same 
topic t^ll. 
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The procedure for performing the above mapping onto an effective Heisenberg Hamil- 
tonian rehes on the constrained density functional theory which allows to obtain the 
ground state energy for a system subject to certain constraints. In the case of magnetic 
interactions, the constraint consists in imposing a given configuration of spin-polarization 
directions, namely, along within the atomic cell i. Note that intracell non-collinearity of 
the spin-polarization is neglected since we are primarily interested in low-energy excitations 
due to intercell non-collinearity. 



Once the exchange parameters Jij are obtained, the spin-dynamics [^|T6|,|17| can be deter- 
mined from the effective Hamiltonian (|Tp and one obtains the result known from spin-wave 
theories of localized ferromagnets: the spin-wave energy E{q) is related to the exchange 
parameters Jij by a simple Fourier transformation 

^(q) = ^ E -^o.-d - exp(zq ■ Ro,)) , (2) 

where Roj=Ro — Rj denote lattice vectors in the real space, q is a vector in the corresponding 
Brillouin zone, and M is the magnetic moment per atom (/i^ is the Bohr magneton). 

There are basically two approaches to calculate the exchange parameters and spin-wave 
energies. The first one which we adopt in the present paper, referred to as the real-space 
approach, consists in calculating directly Jij by employing the change of energy associated 
with a constrained rotation of the spin-polarization axes in cells i and j [0]. In the framework 
of the so-called magnetic force-theorem P,p!8[ the change of the total energy of the system 
can be approximated by the corresponding change of one-particle energies which significantly 
simplifies calculations. The spin- wave energies are then obtained from Eq. (^. In the second 
approach, referred to as the frozen magnon approach, one chooses the constrained spin- 
polarization configuration to be the one of a spin-wave with the wave vector q and computes 
-E'(q) directly by employing the generalized Bloch theorem for a spin-spiral configuration 
||19|| . Like the above one, approach can be implemented with or without using the the 
magnetic force-theorem. Both the real-space approach and the frozen magnon approach 
can be implemented by using either a finite or an infinitesimal rotation, the latter choice is 
usually preferable. The exchange parameters Jij are then obtained by inverting Eq. (^. One 
should also mention a first-principles theory of spin-fluctuations (the so-called disordered 



local moment picture) based on the idea of a generalized Onsager cavity field ^0 . 

The spin- wave stiffness D is given by the curvature of the spin- wave dispersion E{q) at 
q = 0. Although its calculation is in principle straightforward in the real-space approach, 
we shall show that serious difficulties arise due to the Ruderman-Kittel-Kasuya-Yoshida 
(RKKY) character of magnetic interactions in metallic systems. These difficulties have 
been underestimated in a number of previous studies [|],^,|13,14|, and the claimed agreement 



with experiment is thus fortituous. We shall present a procedure allowing to overcome these 
difficulties. In addition, we shall demonstrate that the evaluation of the spin-wave dispersion 
-E'(q) in the real-space approach has to be also done carefully with respect to the convergency 
of results with the number of shells included. 

Finally, to obtain thermodynamic quantities such as the Curie temperature, we apply 
statistical mechanics to the effective Hamiltonian (|l|). In the present paper, we use two 
different approaches to compute the Curie temperature. The first one is the commonly 
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used mean field approximation (MFA). The limitations of this method are well known: it is 
correct only in the limit of high temperatures (above Tq) , and it fails to describe the low- 
temperature collective excitations (spin- waves). The second approach is the Green function 
method within the random phase approximation (RPA) [p| -p6|. The RPA is valid not 
only for high temperatures, but also at low temperatures, and it describes correctly the 
collective excitations (spin- waves). In the intermediate regime (around Tc), it represents a 
rather good approximation which may be viewed as an interpolation between the high and 
low temperature regimes. It usually yields a better estimate of the Curie temperature as 
compared to the MFA. It should be noted, however, that both the MFA and RPA fail to 
describe correctly the critical behavior and yield in particular incorrect critical exponents. 



II. FORMALISM 

The site-off diagonal exchange interactions J^- are calculated using the expression 
1 



Jij = — Ira tiL 
An Jc 



{P}{z)-PKz))gUz){P]{z)-P\{z))g\,{z) dz, (3) 



which is evaluated in the framework of the first-principles tight-binding linear muffin-tin 
orbital method (TB-LMTO) Here tr^, denotes the trace over the angular momentum 
L = (im), energy integration is performed in the upper half of the complex energy plane 
over a contour C starting below the bottom of the valence band and ending at the Fermi 
energy, P[{z) are diagonal matrices of the so-called potential functions of the TB-LMTO 
method for a given spin direction a =t, i with elements P^^l{z) , and g^jiz) are the so-called 
auxiliary Green function matrices with elements glLjL'i^) HH defined as 

la^'i^M^L' = PIl{^) ^L' - S.L.L' . (4) 

We have also introduced the spin-independent screened structure constant matrix Sij with 
elements Snjii which characterizes the underlying lattice within the TB-LMTO approach 

Calculated exchange parameters were further employed to estimate the spin-wave spec- 
trum E[c\) as given by Eq. (^. For cubic systems and in the range of small q we have 

E(q)=Dg^ (5) 

where q = |q|. The spin- wave stiffness coefficient D can be expressed directly in terms of 
the exchange parameters Joj as P] 

where i?oj = |R-oil- The summation in Eq. (|^) runs over all sites but in practice the 
above sum has to be terminated at some maximal value of i?oj = Rmax- There is a lot 
of misunderstanding in the literature as concerns the use of Eq. (B|). Several calculations 



were done with Rmax corresponding to the first few coordination shells p,13,14]. In other 
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calculations the authors realized the problem of the termination of Rmax but they did 

not suggest an appropriate method to perform the sum ((P) in the direct space. We will 
demonstrate that terminating the sum in Eq. (^) after some value of Rmax is fundamentally 
incorrect because it represents a non-converging quantity and we will show how to resolve 
this problem from a numerical point of view. The reason for such behavior is the long-range 
oscillatory character of Jij with the distance Rij in ferromagnetic metals. 

Alternatively, it is possible to evaluate E{q) directly in the reciprocal space ||^ as 

E{ci) = ^ (J(0) - J(q)) , (7) 



dz , 



and to determine the spin-stiffness constant as a second derivative of E{q) with respect to 

q 

Calculated exchange parameters can be also used to determine Curie temperatures of 
considered metals. Within the MFA 

fc.Tr^ = ^E^o. = ^^E^(q), (8) 

where E{q) is the spin-wave energy (|^). We have calculated T^^^^ directly from the expres- 
sion ksT^'^^ = 2Jo/3, where 



J,^J2Jo^ = -^ [ Imtri [(Pj(z) - P^{z)) (gUz) - gUz)) + 

{Pl{z)-P^{z))gUz){Pl{z)-P^{z))gUz)\ dz . (9) 



The expression for the Curie temperature within the GF-RPA approach is [26 



' '"^^E^- (10) 



The integrand in (|TD|) is singular for q = 0. We have therefore calculated Tq^^ using the 
expression 

lim Re (z) , 



The quantity Gm{z) is the magnon Green function corresponding to the dispersion law E{q) 
and it was evaluated for energies z in the complex energy plane and its value for z = was 



obtained using the analytical deconvolution technique |29|. It should be noted that the MFA 
and the RPA differ essentially in the way in which they weight various Jij, namely more 
distant neighbors play a more important role in the RPA as compared to the MFA. It is seen 
from Eqs. (p|,p!OD that T^^"^ and T^^"^ are given as the arithmetic and harmonic averages 
of the spin-wave energies -E(q), respectively, and therefore it holds T^^^ > Tq^^. 
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III. RESULTS AND DISCUSSION 



A. Details of calculations 

Potential functions and Green functions which appear in Eq. were determined within 



assum- 



the non-relativistic TB-LMTO method in the so-called orthogonal representation p8 
ing the experimental lattice constants and the exchange potential in the form suggested by 
Vosko-Wilk-Nusair ||30|. It should be noted that some calculations, in particular for T^^^^, 
were also done using the scalar-relativistic formulation. The contour integral along the path 
C which starts below the lowest occupied band and ends at the Fermi energy (we assume 



zero temperature) was calculated following the scheme described in which employs the 
Gaussian quadrature method. Twenty energy nodes were used on the semi-circle in the 
upper part of the complex energy plane. The integration over the full Brillouin zone was 
performed very carefully to obtain well-converged results even for very distant coordination 
shells (up to 172-nd shell for fee lattice and the 195-th shell for bcc lattice). In particular, we 
have used up to 5 x 10^ k-points in the full Brillouin zone for the energy point on the contour 
C closest to the Fermi energy, and the number of k-points then progressively decreased for 
more distant points, and for points close to the bottom of the band. 



B. Effective Heisenberg exchange parameters 



We will first discuss qualitatively the dependence of Jij on the distance R 



Ri — R,j 



In the limit of large values of Rij the expression (^) can be evaluated analytically by means 



of the stationary phase approximation |3T|]. For simplicity we consider here a single-band 
model but the results can be generalized also to the multiband case (see Ref. |]3^). For a 
large R^j behaves g^^ as 



exp [i(k" ■ R,,- + $")] 



(12) 



where k'^ is the wave vector of energy ii^ in a direction such that the associated group 
velocity Vk-E'^(k) is parallel to Rjj, and denotes a corresponding phase factor. The 
energy integration in (^) yields additional factor of 1/Rij |^] and one obtains 



Jij oc Im 



exp 



^fk' 



kp) ■ Rjj 



(|)T _|_ 



R^ 



(13) 



For a weak ferromagnet both Fermi wave vectors kp and kp are real and one obtains a 
characteristic RKKY-like behavior 



Jij oc 



sm 



(k|. + k^) ■ Ri,- + $T + $i 



p3 



(14) 



i.e., the exchange interaction has an oscillatory character with an envelope decaying as 
1/i?? . On the other hand, for a strong ferromagnet with a fully occupied majority band 
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the corresponding Fermi wave vector is imaginary, namely kp = in^ and one obtains an 
exponentially damped RKKY behavior 

sin(k^ ■ Rjj + + exp(- kI ■ Kjj) 
Jij oc -3 . [Lb) 

The qualitative features of these RKKY-type oscillations of Jij will not be changed in realistic 
ferro magnets. For a weak ferromagnet, like Fe, one expects a pronounced RKKY character 
giving rise to strong Kohn anomalies in the spin-wave spectrum. On the other hand, for 
Co and Ni which are almost strong ferromagnets one expects a less pronounced RKKY 
character, less visible Kohn anomalies in the spin- wave spectrum (see Sec. |III C|) , and faster 
decay of Jij with a distance Rij. It should be noted that due to the sp-d hybridization no 
itinerant ferromagnet is a truly strong ferromagnet. 

The calculated Heisenberg exchange parameters Jij for bcc-Fe, fcc-Co, and fcc-Ni are 
presented in the Table I for the first ten shells. The exchange parameters Jij remain non- 
negligible over a very long range along the [lll]-direction, and change from ferromagnetic 
to antiferromagnetic couplings already for the third nearest-neighbors (NN). In case of Co 
this change appears only for the 4th NN whereas Ni remains ferromagnetic up to the 5th 
NN. It should be noted a short range of Jij for the case of Ni, being essentially a decreasing 
function of the distance with the exception of the second NN. Such behavior is in a qualitative 
agreement with conclusions obtained from the asymptotic behavior of Jij with distance 
discussed above, in particular with the fact that bcc-Fe is a weak ferromagnet while fcc- 
Co and, in particular, fcc-Ni are almost strong ferromagnets. There have been several 
previous calculations of J^'s for Fe and Ni P,|5|,[7|, pJ]j26[| . Present calculations agree well 



with calculations of Refs. [3,^pT| and there is also a reasonable agreement with results of 
Refs. p6|,[7|. It should be mentioned that Jij for both fcc-Co and hcp-Co were determined 
and they agree quite well with each other (see also Table III below). Finally, we have 
also verified numerically the validity of important sum rule, namely Jq = J2i^o Joi- The sum 
fluctuates with the number of shells very weakly for, say, more than 50 shells. 



C. Dispersion relations 

Calculated magnons energy spectra E{q) along the high symmetry directions of the 
Brillouin zone are presented in Figs. |I] a-c together with available experimental data PB|-|55 



We have used all calculated shells to determine E{q), namely 195 and 172 shells for bcc- and 
fcc-metals, respectively. Corresponding plot of E{q) for fcc-Ni exhibits parabolic, almost 
isotropic behavior for long wavelengths and a similar behavior is also found for fcc-Co. On 
the contrary, in bcc-Fe we observe some anisotropy of E{q), i.e., E{q) increases faster along 
the r — N direction and more slowly along the T — P direction. In agreement with Refs. |^,|^ 
we observe a local minima around the point H along T — H and H — N directions in the 
range of short wavelengths. They are indications of the so-called Kohn anomalies which 
are due to long-range interactions mediated by the RKKY interactions similarly like Kohn- 
Migdal anomalies in phonon spectra are due to long-range interactions mediated by Friedel 
oscillations. It should be mentioned that minima in dispersion curve of bcc-Fe appear only if 
the summation in (EI) is done over a sufficiently large number of shells, in the present case for 
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more than 45 shells. A similar observation concerning of the spin-wave spectra of bcc-Fe was 
also done by Wang et al. |26| where authors used the fluctuating band theory method using 
semiempirical approach based on a fitting procedure for parameters of the Hamiltonian. On 
the other hand in a recent paper by Brown et al. above-mentioned Kohn- anomalies in 
the behavior of spin-wave spectra of bcc-Fe were not found, possibly because the spin-wave 
dispersion was obtained as an average over all directions in the q-space. 

Present results for dispersion relations compare well with available experimental data 
of measured spin- wave spectra for Fe and Ni []33|-p5[]. For low-lying part of spectra there 
is also a good agreement of present results for dispersion relations with those of Refs. |^,|^ 
obtained using the frozen magnon approach. There are, however, differences for a higher 
part of spectra, in particular for the magnon bandwidth of bcc-Fe which can be identified 
with the value of E{q) evaluated at the high-symmetry point q = H in the bcc-Brillouin 
zone. The origin of this disagreement is unclear. We have carefully checked the convergence 
of the magnon dispersion laws -E'(q) with the number of shells included in Eq. (0) and it was 
found to be weak for 50 - 70 shells and more. However, if the number of shells is small the 
differences may be pronounced, e.g., our scalar- relativistic calculations give for the bcc-Fe 
magnon bandwiths the values of 441 meV and 550 meV for 15 and 172 shells, respectively. 
The former value agrees incidentally very well with that given in Refs. [^0. On the other 
hand, even small differences in values of E{q) are strongly amplified when one evaluates the 
second derivative of E{q) with respect to q, i.e., the spin-wave stiffness constant. One should 
keep in mind, however, that the above discussion is somehow academic, for it concerns an 
energy region where the adiabatic approximation ceases to be a good one, so that spin-waves 
are non longer well defined because of their strong damping due to Stoner excitations (see 
e.g. [0]. The results of theoretical calculations based upon the adiabatic approximation can 
be thus compared with each other, but not with experimental data. It should be pointed 
out that the influence of deviations in the calculation of magnon spectra for large values of 
q of the Curie temperature is minimized for its RPA value as compared to its MFA value 
(see Eqs. (|,|10D. 



D. Spin- wave stiffness constant 

As was already mentioned the sum in (|^) does not converge due to the characteristic 
RKKY behavior (14) and, therefore, Eq. (^ cannot be used directly to obtain reliable values 



for the spin-wave stiffness constant. This is demonstrated in Fig. |^ where the dependence 
of calculated spin-wave stiffness constants on the parameter Rmax in Eq- (ID is plotted. The 
oscillatory character of D versus Rmax persists for large values of Rmax for the case of bcc-Fe 
and even negative values of spin-wave stiffness constants were obtained for some values of 
Rmax- To resolve this difficulty we suggest to regularize the expression (^ by substituting 
it by the formally equivalent expression which is, however, numerically convergent 

D = lim D{7]) , 



D{rj) = \im ^ ^ Joj RIj exp{-r]Roj/a) . (16) 
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The quantity rj plays a role of a damping parameter which makes the sum over Rij absolutely 
convergent as it is seen from Fig. |^. The quantity 0(7]) is thus an analytical function of the 
variable t] for any value 77 > and can be extrapolated to the value 77 = 0. To show that the 
limit for 77 — >■ is indeed finite and that our scheme is mathematically sound, let us consider 
as an example a typical RKKY interaction J{R) oc sm{kR + For large R we can 

employ Eq. (|l^ and substitute the sum in @ by a corresponding integral. We obtain 



lim DM oc 47r H ^^'^(^-^ + dR = - cos(<l>) sUkRo) - sin($) ci(kRo) , (17) 

'7^0 jRg Ri 



where si and ci are integral sine and cosine, respectively. The integral is indeed finite. 

We therefore perform calculations for a set of values rj G {rjminyVmax) for which D{ri) is 
a smooth function with a well pronounced limit for large Rmax- The limit rj = is then 
determined at the end of calculations by a quadratic least-square extrapolation method. 
Typically, 5-15 values of rj was used for rjmin ~ 0.5 — 0.6 and rjmax ~ 0.9 — 1.2 with a relative 
error of order of a few per cent. In calculations we have used i?max=7a for fee and 9a for bcc, 
where a denotes the corresponding lattice constant. It should be noted a proper order of 
limits in Eq. (|16]), namely first evaluate a sum for large Rmax and then limit r] to zero. The 
procedure is illustrated in Fig. ^. The results for spin stiffness coefficient D calculated in this 



way are summarized in Table II together with available experimental data |36-3S]. There is 
a reasonable agreement between theory and experiment for bcc-Fe and fcc-Co but the values 
of spin-wave stiffness constant are considerably overestimated for fcc-Ni. It should be noted 
that measurements refer to the hcp-Co while the present calculations were performed for 
fcc-Co. A similar accuracy between calculated and measured spin-wave stiffness constants 
was obtained by Halilov et al. @] using the frozen-magnon approach. Our results are also 
in a good agreement with those obtained by van Schilfgaarde and Antropov ||^ using the 
spin-spiral calculations to overcome the problem of evaluation of D from Eq. (j^). On the 
other hand, this problem was overlooked in Refs. P jr3|JT^ so that a good agreement of D, 



calculated for a small number of coordination shells, with experimental data seems to be 
fortituous. Finally, results of Brown et al. |^ obtained by the layer Korringa-Kohn-Rostoker 
(KKR) method in the frozen potential approximation are underestimated for all metals and 
the best agreement is obtained for Ni. 



E. Curie temperature 



Several attempts have been made to evaluate Curie temperatures of magnetic transition 
metals most of them based on the MFA. The MFA as a rule overestimates 

values of Curie temperatures (with exception of fcc-Ni with values substantially underesti- 
mated). We will show that an alternative method based on the Green function approach 



in the framework of the RPA p4-p5| can give a better agreement with experimental data. 
The RPA Curie temperatures were calculated from Eq. (0) by employing the method of 
analytical deconvolution [^]. In order to test the accuracy of this procedure we compare 
the present numerical results for the ratio j"^^^^ /T^^^ obtained for the nearest-neighbor 
Heisenberg model with the exact results p2|j25[] : we obtain 1.33 (fee) and 1.37 (bcc) as 
compared to exact values 1.34 and 1.39, respectively, i.e., a numerical procedure agrees with 
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exact results within one per cent accuracy. Calculated values of Curie temperatures for both 
the MFA and RPA as well as corresponding experimental data are summarized in Table II. 
Mean-field values of Curie temperatures are overestimated for Fe and Co, but underesti- 
mated for Ni in agreement with other calculations @J^. On the other hand, the results 
obtained using the RPA approach are in a good agreement with experiment for both fcc-Co 
and bcc-Fe, while the results for fcc-Ni are even more underestimated. This is in agreement 
with the fact mentioned in Sec. |T|, namely that Tq^^ < T^^^. The present results for Fe 
and Ni are in a good agreement with results of Ref. using the spin-fluctuation theory 
and an improved statistical treatment in the framework of the Onsager cavity-field method. 

The calculated ratio J'A'/fa^j'Rpa jg 3^ 49^ 1.25, 1.13 for bcc-Fe, fcc-Co, and fcc-Ni, respec- 
tively. The values differ from those obtained for the first-nearest neighbor Heisenberg model 
due to non-negligible next-nearest neighbors in realistic ferromagnets and their oscillatory 
behavior with the shell number. 

The last point concerns the relevance of relativistic corrections for the evaluation of the 
exchange parameters and related quantities. The simplest quantity to evaluate is the MFA 
value of the Curie temperature (see Eq. (||)). Results for ferromagnetic metals (including 
hcp-Co) are summarized in Table III by comparing the non-relativistic and scalar-relativistic 
values. One can conclude that scalar-relativistic corrections are not important for fcc-Co 
and hcp-Co but their effect is non-negligible for fcc-Ni and bcc-Fe. The scalar-relativistic 
corrections generally shifts sp-bands downwards as compared to the d-band complex while 
the changes of magnetic moments are generally very small (a similar exchange splitting). 
One can thus ascribe above changes mostly to the modifications of the density of states at 
the Fermi energy (the site-diagonal blocks of the Green function in Eq. (^). Results also 
show only a weak dependence of the calculated T^^^ on the structure (hcp-Co vs fcc-Co). 



F. Comparison between the real-space and frozen magnon approaches 

The real-space and frozen magnon approaches are formally equivalent to each other. The 
quantities that are directly calculated (the Jj^s in the former case, the -E'(q)'s in the latter) 
are related to each other by a Fourier transformation. Therefore, the pros and cons of both 
approaches concern mainly their computational efficiency. 

The computational effort needed to obtain one Jij parameter within the real-space ap- 
proach is approximately the same as to compute one magnon energy -E(q) within the frozen 
magnon approach: in both fine Brillouin zone integration is required. 

Therefore, it is quite clear that if one is primarily interested in spin-wave dispersion 
curves (for a moderate number of q points), or in the spin- wave stiffness D, the frozen 
magnon approach is superior, for it does require to perform a Fourier transformation and 



the delicate analysis explained in Sec. |IIID| . We have shown, however, although less direct 
and computationally more demanding, the real-space approach performs well also. 

On the other, if one is interested in the Curie temperature, the real-space approach is 
more efficient. This obvious if ones uses the mean- field approximation. Indeed, T^'^^ is 
obtained from a single real-space calculation, by using the sum rule (P), whereas many 
i?(q)'s are needed to obtain T^^^^ from Eq. (^ within the frozen magnon approach. Also 
if one uses the RPA, the real-space approach is more efficient. For both approaches, the 
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integral over q in Eq. (|TOD needs to be performed accurately, with paying great attention to 
the divergence of the integrant at q = 0. A very high density of q points is required there, in 
order to have a satisfactory convergence. Within the frozen-magnon approach, each of the 
-E'(q)'s requires the same computational effort. In contrast, within the real-space approach, 
less than 200 Jj/s are sufficient to obtain a parametrization of E{q) over the full Brillouin 
zone, which considerably reduces the computational effort. 

A further very important advantage of the real-space approach is its straightforward 
application to systems with broken translational symmetry like, e.g., substitutional alloys, 
surfaces, overlayers, and multilayers. This is an important advantage keeping in mind the 
relevance and yet not fully understood character of exchange interactions at metal interfaces 
and surfaces. The reciprocal-space approach can be applied to ideal surfaces, but it is 
numericalUy demanding, and its application to system with substitutional disorder and/or 
to finite magnetic clusters is practically impossible. Finally, the dependence of exchange 
parameters Jij on the distance also gives an important information about the nature of 
the magnetic state (RKKY-like interactions) and this dependence is again straigtforwardly 
determined by the real-space method while in the reciproacal-space method Ji/s have to be 
determined by inverting Eq. (0). 

IV. SUMMARY 

We have calculated Heisenberg exchange parameters of bcc-Fe, fcc-Co, and fcc-Ni in real 
space from first-principles by employing the magnetic force theorem. We have determined 
dispersion curves of magnetic excitations along high-symmetry directions in the Brillouin 
zone, spin-wave stiffness constants, and Curie temperatures of considered metals on the 
same footing, namely all based on calculated values of exchange parameters Jij. Dispersion 
curves of bcc Fe exhibit an anisotropic behavior in the range of long wavelengths, with 
peculiar minima for short wavelengths in the [100]-direction which are due to a relatively 
strong exchange oscillations in this metal. We have presented a method of evaluation of the 
spin-wave stiffness constants which yields converged values, in contrast to previous results in 
the literature. Calculated spin-wave stiffness constants agree reasonably well with available 
experimental data for Co and Fe, while agreement is rather poor for Ni. Present calcula- 
tions agree also well with available experimental data for magnon dispersion law of bcc-Fe. 
We have also evaluated Curie temperatures of metals in question using the mean-field ap- 
proximation and the Green function random phase approximation. We have found that in 
the latter case a good agreement with the experiment is obtained for Co and Fe, while less 
satisfactory results are obtained for Ni, where the role of the Stoner excitations is much 
more important as compared to Co and Fe. In addition, the adiabatic approximation is less 
justified for Ni, and, possibly, correlation effects beyond the local density approximation 
play the more important role for this ferromagnet. 

In conclusion, we have demonstrated that the real-space appraoch is able to determine the 
low-lying excitations in ferromagnetic metals with an accuracy comparable to the reciprocal- 
space approach. This justifies the use of the real-space approach for more interesting and 
complex systems with violated translational symmetry where the reciprocal-space approach 
is of the limited use. In particular, a first promising application of the real-space approach 
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to the problem of the oscillatory Curie temperature of two-dimensional ferromagnets has 
been recently published [^|. 
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TABLES 



TABLE L Effective Heisenberg exchange parameters Jqj for ferromagnetic Fe, Co, and Ni 
for the first 10 shehs. Quantities Rqj and Ns denote, respectively, shell coordinates in units of 
corresponding lattice constants and the number of equivalent sites in the shell. 





Fe (bcc) 






Co (fee) 






Ni (fee) 




Roj 




Joj [mRy] 






Joj [mRy] 






Joj [mRy] 


fill) 


8 


1.432 




12 


1.085 




12 


0.206 


(100) 


6 


0.815 


(100) 


6 


0.110 


(100) 


6 


0.006 


(110) 


12 


-0.016 




24 


0.116 




24 


0.026 


(111) 

\2 2 2l 


24 


-0.126 


(110) 


12 


-0.090 


(110) 


12 


0.012 


(111) 


8 


-0.146 


(i|o) 


24 


0.026 


(i|o) 


24 


0.003 


(200) 


6 


0.062 


(111) 


8 


0.043 


(111) 


8 


-0.003 


f331x 
\2 2 2> 


24 


0.001 




48 


-0.024 




48 


0.007 


(210) 


24 


0.015 


(200) 


6 


0.012 


(200) 


6 


-0.001 


(211) 


24 


-0.032 


(ifo) 


12 


0.026 


(ifo) 


12 


-0.011 


('3 3 3\ 
^2 2 21 


8 


0.187 


(2H) 


2 4 


0.006 


(2^) 


24 


0.001 



TABLE IL Calculated spin-wave stiffness constants {Pth) and Curie temperatures {Tq and 
Tq^^) and their comparison with experimental values D^x and T^^ . 



Metal 


Ah[meV-A'] 


Z)ex[meV-A'] 




r^^^[K] 


wm 


Fe (bcc) 


250 ±7 


280^330'^ 


1414 


950 ±2 


1044 - 1045 


Co (fee) 


663 ±6 


580"'^510'^ 


1645 


1311 ±4 


1388 - 1398" 


Ni (fee) 


756 ± 29 


555^422^ 


397 


350 ±2 


624 - 631 



"Data refer to hep Co at 4.2 K. 

''Neutron scattering measurement at 4.2 K [p^] . 

^Magnetization measurement [^ at 4.2 K. 

'^Neutron scattering measurement extrapolated to K 
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TABLE III. Calculated Curie temperatures of ferromagnetic metals in the mean-field approx- 
imation for non-relativistic (nr) and scalar-relativistic (sr) cases. 



Metal 




rj-[K] 


Fe (bcc) 


1414 


1335 


Co (fee) 


1645 


1651 


Co (hep) 


1679 


1673 


Ni (fee) 


397 


428 
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FIGURES 



FIG. 1. Magnon dispersion relations along high-symmetry lines in the Brillouin zone: (a) bcc-Fe 
(experiment: Ref. [33], 10 K, filled circles and Ref. [35], Fe(12 % Si), room temperature, empty 
squares); (b) fcc-Co; and (c) fcc-Ni (experiment: Ref. [34], room temperature, empty circles) . Lines 
are calculated results. 

FIG. 2. Spin-wave stiffness constants calculated from Eq. ( p!6D as a function of Rmax (in units 
of lattice constants) for fcc-Ni (full line), fcc-Co (short dashes), and bcc-Fe (long dashes). 

FIG. 3. Spin-wave stiffness of fcc-Ni calculated from Eq. (^) as a function of Rmax (in units of 
lattice constant) for various values of the damping factor ry. 

FIG. 4. Spin-wave stiffness coefficients D['q) for bcc-Fe (empty squares), fcc-Co (empty trian- 
gles), and fcc-Ni (empty circles) as a function of the parameter r/ and extrapolated values for r/ = 
(filled symbols). The solid line indicates the quadratic fit function used for extrapolation. 
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